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ABSTRACT 

Molecular  beam  techniques  and  a  time-af-flight  velocity  distribution  detector 
system  were  used  to  investigate  the  effects  of  condensation  in  free-jet  expansions  of 
argon.  A  l-in.-diam  skimmer  cooled  to  22° K  was  used  to  eliminate  skimmer  interaction. 
The  data  reported  are  for  170,  300,  and  600°  K  argon  expanded  into  a  vacuum  from 
source  pressures  of  from  10  to  10,000  torr  and  with  a  nozzle-skimmer  separation  distance 
of  650  nozzle  diameters.  The  data  revealed  that  the  mean  velocity  of  the  monomer 
species  can  be  greater  than  predicted  by  cpT0  =  V4mv^  when  condensation  occurs.  This  is 
hypothesized  to  be  the  result  of  the  heat  of  condensation  being  added  to  the  flow  field. 
Also,  when  condensation  occurs,  the  static  temperature  increases.  Terminal  speed  ratios  as 
high  as  27.5  were  measured  during  these  tests. 
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NOMENCLATURE 

b0  Constant,  b0  =  2/3  ir  N0d3 

C  Constant 

Cp  Specific  heat  at  constant  pressure 

D„  Orifice  diameter 

d  Molecular  diameter 

Fr  Free  energy  of  a  drop  of  radius,  r 

F«.  Free  energy  of  a  large  liquid  surface 

fr  Fugacity  of  a  drop  of  radius,  r,  at  gas  and  liquid  phase  interface 

f„  Fugacity  of  a  large  liquid  surface  at  gas  and  liquid  phase  interface 

K  Initial  relative  kinetic  energy  or  the  total  energy  of  a  center  of  mass  system 

Kns  Knudsen  number  evaluated  at  the  source 
k  Boltzmann's  constant 

L  Flight  length  from  chopper  to  detector 

M  Mach  number 

MK  Condensation  Mach  number 

Mj  Terminal  Mach  number 

m  Mass  of  a  molecule 

N  Total  beam  flux 

N0  Avagodro’s  number 

P  Pressure 

PB  Sum  of  pressures  produced  by  beam  and  background  gas 

PBC  Pressure  of  background  gas  in  chamber 

pK  Actual  pressure  of  condensation 
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p0  Source  pressure 

Py  Vapor  pressure  of  a  drop  of  radius,  r 

P«,  Vapor  pressure  of  a  flat  liquid  surface 

R  Gas  constant 

r  Molecular  radius  or  molecular  separation  in  Eqs.  (8)  and  (9) 

t*  Critical  radius  of  a  drop  which  must  be  present  for  condensation  to  occur 

S<j  (t)  Signal  distribution  function 

SR  Speed  ratio 

T  Temperature 

Tic  Actual  temperature  of  condensation 
T0  Source  temperature 

Ts  Static  temperature 

t  Time 

tc  Time  delay  between  photoelectric  trigger  and  beam  molecules  leaving  chopper 

tp  Waveform  eductor  time  delay 

tm  Time  of  maximum  signal 

t„  Time  associated  with  flow  velocity,  U0 

t,  Spectrometer  time  delay— the  time  between  beam  ionization  and  detection 

U  Centrifugal  potential 

U0  Flow  velocity  of  the  beam 

v  Thermal  velocity 

vra  Velocity  corresponding  to  the  time  of  maximum  signal 

vm  p  Most  probable  velocity 

W  Thermodynamic  probability 
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X2B  Concentration  of  bound  dimers 

a  Lower  limit  of  molecular  separation  for  which  repulsive  potential,  0(r),  is  2ero 

ai  Beam  flux  gage  calibration  constant 

0  Constant  0  =  m/2kT 

y  Specific  heat  ratio 

At  “On”  time  of  a  single  waveform  eductor  module 

6  Total  depth  of  potential  well 

Xo  Mean  free  path  at  source  conditions 

r}  Surface  tension 

p  Density 

Pl  Density  of  liquid  phase 

p0  Density  at  source  conditions 

0  Gas  supersaturation 

0k  Gas  supersaturation  at  which  condensation  occurs 

0(r)  Interaction  potential 
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SECTION  I 
INTRODUCTION 

In  the  normal  theory  of  condensation  in  large-scale  jets,  the  condensation  line 
uniquely  defines  the  pressure  and  temperature  conditions  which  are  necessary  before  a 
pure,  neutral  gas  can  condense.  The  normal  theory  is  adequate  if  a  small  mean-free  path 
exists  between  particles  in  an  expanding  jet  since  the  interaction  between  particles  results 
in  a  quasi-equilibrium  condition. 

With  the  advent  of  the  space  age,  it  has  become  apparent  that  the  normal  theory 
is  inadequate  for  analyzing  condensation  effects  in  a  free-jet  expansion  (i.e.,  expansion 
from  high-pressure  reservoir  conditions  into  a  near  vacuum).  The  greater  mean-free  path 
between  particles  results  in  a  much  greater  time  rate  associated  with  the  condensation 
process,  and  the  quasi-equilibrium  assumptions  are  no  longer  valid. 

A  better  understanding  of  condensation  and  molecular  cluster  formation  in  a 
free-jet  expansion  is  needed.  These  phenomena  are  important  considerations  in 
space-related  problems,  such  as  determining 

1 .  The  gas  dynamic  characteristics  of  highly  expanded  rocket  gas  plumes 
from  space  propulsion  systems, 

2.  The  effect  of  the  forces  generated  by  molecular  clusters  when  these 
clusters  impinge  on  a  body  immersed  in  a  rocket  plume,  and 

3.  The  expansion  characteristics  of  reactive  gases  such  as  nitrous  oxide 
(NO)  which  are  used  in  high  altitude  research. 

An  approach  to  analyzing  the  condensation  phenomena  in  a  free-jet  expansion  is 
to  allow  the  expansion  to  proceed  to  its  terminal  Mach  number,  isolate  the  core  of  the 
expanded  jet,  and  measure  the  velocity  distribution  of  the  different  size  molecular 
clusters  in  this  core.  The  measurement  of  the  relative  abundances  and  velocity 
distributions  of  molecular  dusters  in  the  condensation  region  provides  additional  data 
which  can  be  used  in  the  analysis  of  the  mechanism  of  condensation. 

Concurrent  developments  in  aerodynamic  molecular  beam  technology  and  the 
theory  of  condensation  in  a  free  jet  during  the  last  two  decades  allow  the  approach 
outlined  above  to  be  used.  These  developments  fall  into  three  broad  categories: 

1.  Improvements  in  distinguishing  between  the  effects  of  the  molecular 
beam  and  the  background  gas  present  in  the  molecular  beam  chamber, 

2.  The  development  of  methods  for  analyzing  time-of- flight  (TOF)  data 
and  of  obtaining  simple  moments  (i.e.,  weighted  integrals)  based  on 
information  which  is  directly  available  in  terms  of  TOF 
measurements,  and 
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3.  The  development  of  a  qualitatively  valid  model  of  the  condensation 
process  in  a  free  jet. 

Concurrent  with  the  improvements  in  aerodynamic  molecular  beam  technology, 
various  Investigators  contributed  bits  and  pieces  of  evidence  which  indicate  that  the 
formation  of  molecular  clusters  is  an  intrinsic  part  of  the  condensation  process.  The 
development  of  this  theory  is  chronologically  developed  in  the  following  paragraphs. 

If  a  gas  in  a  cloud  chamber  apparatus  is  expanded  into  a  region  of  lower  pressure, 
the  degree  of  expansion  required  to  achieve  condensation  is  determined  by  the  point  of 
intersection  of  the  pressure- temperature  adiabat  and  the  condensation  line.  The 
condensation  line  uniquely  defines  the  pressure  and  temperature  conditions  under  which 
a  pure  neutral  gas  will  condense.  Normally  the  gas  does  not  condense  when  it  reaches  the 
saturated  state,  but  instead  it  requires  a  certain  degree  of  supersaturation  to  bring  about 
droplet  formation.  If  Pk  and  Tg  are  the  actual  pressure  and  temperature  of 
condensation,  the  appropriate  supersaturation,  0k  >  is  defined  by  Pk/p«  (Ref.  I). 

The  relationship  between  the  fugacity  of  a  drop  in  equilibrium  with  a  gas  over  a 
spherical  drop  surface  of  radius,  r,  and  the  fugacity  of  a  large  liquid  surface  in 
equilibrium  with  the  same  gas  can  be  expressed  by 

Fr  -  Fw  «  RT  -Ji.  =  (1) 

feo  r 


as  developed  in  Ref.  2. 

For  a  perface  gas,  the  ratio  of  vapor  pressures  and  the  fugacity  ratio  are  identical. 
By  substituting  vapor  pressures  for  fugacities  and  manipulating  the  equation,  the 
Kelvin-Helmholtz  equation  is  produced: 


P„  r  Pl>T 

where  ij  is  the  surface  tension,  k  is  the  Boltzmann  constant,  Pl  is  the  liquid  density,  and 
m  is  the  mass  of  the  molecule. 

Since  the  ratio,  pt/pm,  is  a  measure  of  the  supersaturation,  0,  it  follows  from  Eq. 
(2)  that  for  every  value  of  0,  there  exists  some  radius,  r,  of  liquid  drops  which  would  be 
just  the  right  size  to  remain  in  equilibrium  with  the  vapor. 


When  condensation  actually  does  occur,  then 


„  pr*  .  2  71  m 

-  =  En.  <f>\  -  - 


pi_kT 


(3) 


The  term,  r*,  is  defined  as  the  “critical  radius.”  This  radius  represents  the  minimum  size 
of  drop  which  must  be  present  for  condensation  to  occur. 
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The  concept  of  a  “critical  radius”  leads  to  the  conclusion  that  a  supersaturated 
vapor  cannot  condense  until  it  contains  certain  “condensation  nucleii.”  In  1 897,  Wilson 
suggested  that  a  pure  neutral  gas  at  a  sufficiently  large  value  of  0  would  contain  these 
nucleii  in  the  form  of  small  molecular  clusters.  Wilson  postulated  that  intermolecular 
forces  gave  each  molecular  cluster  a  finite  lifetime  during  which  it  could  either  grow  by 
assimilation  of  further  molecules  or  decay  through  evaporation.  The  size  of  the  molecular 
cluster  would  determine  which  of  the  two  processes  was  more  probable.  He  further 
postulated  that  the  probabilities  of  growth  or  decay  must  be  equal  for  an  assembly  of 
critical  size  so  that  stable  equilibrium  with  the  vapor  is  established.  This  growth  to  a 
critical  radius  has  come  to  be  called  “self-nucleation”  and  the  sub-critical  and  critical 
assemblies  are  denoted  as  “condensation  embryos”  and  “condensation  nucleii,” 
respectively. 

The  early  cloud  chamber  research  introduced  the  basic  concept  that  for 
condensation  to  occur,  an  assembly  of  “critical  radius”  condensation  embryos  must  be 
present.  This  idea  proved  valid  in  condensation  studies  of  expanding  jets  under 
continuum  flow  conditions. 

In  a  free-jet  expansion,  the  greater  mean-free  path  between  particles  results  in  a 
much  greater  time  rate  associated  with  the  condensation  process,  and  the 
quasi-equilibrium  assumptions  are  no  longer  valid.  The  physics  of  the  problem,  however, 
would  still  lead  to  the  belief  that  “critical  radius”  molecular  clusters  must  exist 
somewhere  in  the  free-jet  flow  field  before  droplet  formation  can  occur. 


Leckenby,  et  al.  (Ref.  1 ),  proposed  that  the  formation  of  molecular  clusters  in  a 
free  jet  was  caused  by  a  condensation  process  brought  about  by  adiabatic  cooling  as  the 
gas  expanded  from  the  source  orifice.  In  the  continuum  region  of  the  expansion,  the  gas 
gains  linear  kinetic  energy  at  the  expense  of  its  enthalpy  and  accelerates  to  supersonic 
Mach  numbers  while  the  temperature,  pressure,  and  density  all  decrease.  The  relationship 
between  these  parameters  is  given  by  the  well-known  equation: 


—  =  11  +  -(y  -  l)M2r* 

T„  2  ■ 


(y-U/y 


Leckenby  expressed  the  line  of  a  saturated  vapor  pressure  curve  as 


tn  -  A  -  B/T 


(4) 


(5) 


where  p_  is  the  vapor  pressure  over  a  flat  surface,  and  A  and  B  are  constants  for  a 
particular  vapor.  The  vapor  does  not  condense  on  reaching  the  saturated  state  but 
requires  a  certain  degree  of  supersaturation  before  droplet  formation  occurs.  The 
appropriate  supersaturation  was  previously  defined  as 


(6) 
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Leckenby,  therefore,  expressed  the  equation  of  the  condensation  line  by  manipulating  the 
equation  for  the  saturated  vapor  line,  Eq.  (5),  as  follows: 


Equation  (4)  could  now  be  used  to  determine  the  Mach  number  at  the  point  of 
condensation : 

— y  ~j 

in  h  +  —  (y  -  l)Mk2  y~l  =  A  -  —  [1  +  —  (y  -  1)  Mk2]  (8) 
L<^k  '  2  -I  Tp  2 

Equation  (8)  implies  that  early  condensation  (i.e.,  small  Mr)  can  be  affected  by 
increasing  p0 ,  decreasing  Tc ,  and  using  a  gas  with  a  large  7  value. 

Leckenby  realized  that  a  free-jet  expansion  involves  a  transition  from  continuum 
flow  to  free  molecular  flow  during  which  the  adiabatic  cooling  ceases,  and  after  which 
temperature  and  flow  velocity  remain  constant.  The  final  temperature  and  flow  velocity 
determine  the  terminal  Mach  number,  My .  He,  therefore,  established  the  condensation 
criteria  that  MT  must  be  greater  than  MK  before  condensation  conditions  are  possible. 

Based  on  the  relationships  developed  in  Appendix  II,  Leckenby  concluded  that 
the  terminal  Mach  number  depends  on  the  reservoir  conditions  (pD,  T0)  and  7  in  a 
directly  contrary  sense  to  the  condensation  Mach  number,  Mr.  He  concludes  (Ref.  1) 
that  “high  reservoir  pressures,  low  initial  temperatures,  and  the  use  of  a  gas  with  a  high  7 
are  conducive  to  well-advanced  expansions  (large  My)  and  early  condensation  (low  Mr).” 

Lenkenby  also  noted  (Ref.  1)  that  the  observation  of  molecular  clusters  is  so 
closely  tied  to  the  expected  occurrence  of  condensation  in  the  gas  “that  it  is  difficult  to 
avoid  the  inference  that  the  one  is  a  direct  manifestation  of  the  other.  The  results  show 
that,  in  general,  the  clusters  appear  in  greater  abundance  and  at  lower  pressures  for  gases 
which  have  high  critical  temperatures  (i.e.,  are  easily  liquified  by  cooling),  except  insofar 
as  the  expansion  is  modified  by  the  variable  7.” 

The  existence  of  dimers  (two  atoms  bound  together  by  van  der  Waals  interaction) 
is  well  documented  in  the  literature  (Refs.  1,  3,  4,  and  5). 

Jeans  (1904)  envisioned  the  formation  of  dimers  as  resulting  from  three-body 
collisions  in  which  the  act  of  binding  occurs  when  excess  energy  is  removed  by  some 
third  body. 

Leckenby  and  Robbins  (Ref.  1)  proposed  in  1965  that  dimers  are  the  smallest 
condensation  embryos  present  in  a  gas.  These  dimers  exist  in  all  gases  if  the  temperature 
is  low  enough  for  the  attractive  van  der  Waals  interaction  to  be  significant  and  when  the 
pressure  is  high  enough  to  produce  a  significant  number  of  three-body  collisions.  The 
dimers  exist  independent  of  whether  the  gas  is  in  a  saturated  or  unsaturated  condition. 


4 


AEDC-TR-70-142 


Stogryn  and  Hirschfelder  (Ref.  6)  have  shown  that,  at  temperatures  less  than  the 
critical  temperature,  the  dominant  contribution  to  the  second  virial  coefficient  is  from 
bound  dimers.  They  have  computed  that  two  types  of  dimers  can  exist:  (1)  bound  dimers 
and  (2)  metastable  dimers. 

If  an  interaction  potential,  $(r),  is  considered,  it  may  be  shown  that,  by  applying 
Maxwell-Boltzmann  statistics  to  the  system  of  interacting  atom  pairs,  the  probability  of 
finding  a  pair  with  separation  between  r  and  (r  +  dr)  and  with  relative  kinetic  energy 
between  U  and  (U  +  dU)  is: 


d2W  =  A  exp  ^  7  -rl\J  U5*  dU  r2  dr  (9) 

where  A  is  normalizing  constant  (Ref.  2). 

Since  the  long  range  forces  between  molecules  are  attractive  (Ref.  7),  it  follows 
that  0(r)  is  negative  over  a  wide  range  of  molecular  separations.  Hirschfelder,  et  al.  (Ref. 
8),  concluded  that,  for  certain  pairs  of  atoms,  the  total  energy  U  +  tf>(r)  would  be 
negative.  This  combination  was  defined  as  “bound”  since,  in  the  absence  of  external 
interactions,  each  atom  is  trapped  in  the  potential  well  of  its  partner.  The  result  of  this 
type  of  bond  is  that  the  two  atoms  would  orbit  indefinitely  about  their  common  center 
of  mass  if  no  collisions  occurred.  In  reality,  the  effect  of  collisions  would  soon> result  in 
the  pair  either:  (1)  absorbing  enough  energy  from  collisions  to  separate,  or  (2)  capturing 
the  colliding  atom  and  forming  a  three-particle  cluster.  The  occurrence  of  either  event 
would  be  determined  by  the  relative  positions  of  the  bound  pair  at  the  time  of  collision. 

Stogryn,  et  al.  (Ref.  6),  calculated  the  concentration  of  bound  dimers,  X2B,  for  a 
closed  system  in  equilibrium  based  on  the  Lennard  Jones  potential: 

-4((t)'!-(t)1  (10) 

where  e  is  the  maximum  depth  of  the  potential  well  and  a  is  the  lower  limit  of  r  for 
which  4>(t)  is  zero.  They  showed  that 


where  „ 

b0  -  -3-  »  N„d3 

v0  is  the  molar  volume  at  standard  conditions,  and  p  is  the  pressure  in  atmospheres. 
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The  essential  features  of  Eq.  (11)  are  that  the  concentration  of  bound  dimers, 
Xjb  >  varied  linearly  with  pressure  and  approximately  as  T"s/J.  It  can  also  be  concluded 
that  gases  with  large  attractive  potentials  (large  e)  most  readily  form  dimers. 

Reference  6  notes  that  in  addition  to  bound  dimers  which  must  meet  the 
criterion  of  (U  +  0(r))  less  than  zero,  another  situation  exists  which  can  result  in  the 
formation  of  “metastable  bound”  dimers. 


If  0(r)  varies  more  rapidly  than  1/r2,  orbiting  pairs  may  exist  even  if  the  total 
energy  is  positive.  This  situation  results  because  the  radial  action  of  the  interacting 
particle,  when  reduced  to  a  one-dimensional  problem,  is  governed  by  both  the  actual 
potential,  0(r),  and  the  “centrifugal  potential”  which  may  be  written  Kb2/r2  where  K  is 
the  initial  relative  kinetic  energy  and  b  is  the  impact  parameter.  The  combination  of  0(r) 
and  the  centrifugal  potential  (which  is  always  positive)  has  a  positive  maximum 
somewhere  in  the  attractive  region  of  0(r).  The  height  and  position  of  this  maximum  are 
determined  by  the  parameters  K  and  b  as  shown  in  Fig.  1  (Appendix  I).  A  potential 
“ditch”  exists  between  this  “hump”  and  the  repulsive  part  of  0(r)  which  confines  the 
radial  oscillations  of  any  atom  pair  whose  total  energy  is  less  than  the  hump  height. 
■These  pairs  fall  into  two  classes,  depending  on  whether  their  total  energy  is  positive  or 
negative. 


and 


Classical  mechanics  makes  no  distinction  between  pairs  where 
U  +  d(r)  <  0  (Bound  dimer  criterion) 

U  +  d(r)  >  0  (Metastable  dimer  condition) 


Leckenby,  et  al.  (Ref.  1),  points  out  that  from  the  standpoint  of  wave  mechanics,  “pairs 
iwith  positive  energy  (metastable  dimers)  are  liable  to  disassociate  by  tunnel  effect  leakage 
through  the  hump.” 

<  Reference  6  gives  a  detailed  treatment  of  bound  and  metastably  bound  dimers  in 

;gases.  These  results  show  that  the  fraction  of  metastably  bound  dimers  is  normally 
somewhat  less  than  one-third  of  the  total  dimer  abundance. 


t  The  central  point  to  be  made  in  this  rather  lengthy  development  is  that  theory 
.leads  to  the  recognition  of  two  categories  of  dimers  which  both  contribute  significant 
^amounts  to  the  total  dimer  abundance.  Classical  statistical  mechanics  cannot  distinguish 
between  bound  and  metastably  bound  dimers.  The  metastably  bound  dimer,  which  has  a 
: positive  total  energy,  would  be  assumed  to  have  a  greater  tendency  to  disassociate  upon 
^collision,  and  it  would  also  tend  to  have  a  higher  average  energy  content. 

In  summary,  aerodynamic  molecular  beam  technology  has  progressed  to  the  point 
where  it  can  provide  useful  data  about  the  condensation  process  in  a  free  jet.  The  degree 
of  correlation  between  these  data  and  the  existing  free-jet  condensation  theory  provides 
the  ‘‘building  blocks”  which  are  necessary  for  a  more  complete  understanding  of  the 
condensation  process. 
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This  investigation  will  experimentally  determine  the  velocity  distribution  and 
relative  abundances  of  the  different  mass  number  argon  clusters  in  the  condensation 
region  of  a  free-jet  expansion  over  a  source  pressure  range  of  from  10  to  10,000  torr  and 
source  temperatures  of  170,  300,  and  600°  K. 

SECTION  II 
APPARATUS 

Z1  AERODYNAMIC  MOLECULAR  BEAM  SYSTEM 

The  AEDC  aerodynamic  molecular  beam  system  (Fig.  2)  is  a  stainless  steel 
cylinder,  3  ft  in  diameter  by  6-1/2  ft  long,  which  is  divided  into  three  sections  by  two 
removable  bulkheads.  Vacuum  conditions  are  produced  and  maintained  in  the  cell  by  oil 
diffusion  pumps,  20° K  gaseous-helium-cooled  cryoliners,  and  77° K  liquid-nitrogen-cooled 
cryoliners.  The  total  pumping  speed  for  air  is  in  excess  of  500,000  liters/sec.  This  system 
maintains  pressures  of  1  x  Iff4 ,  1  x  Iff7,  and  3  x  Iff8  torr  in  the  nozzle  discharge, 
collimation,  and  test  sections  of  the  chamber,  respectively. 

The  apparatus  used  for  generating  the  molecular  beam  during  this  series  of  tests 
consisted  of  a  4-mm-diam  collimating  orifice,  a  2.54-cm-diam  cryogenically  cooled 
skimmer,  a  collimation  distance  of  44  cm  a  0.34-mm-diam  nozzle,  and  a  nozzle-skimmer 
separation  distance  of  22  cm.  A  complete  description  of  the  beam  system  and  its 
performance  is  given  in  Ref.  9. 

2.2  BEAM  SOURCE 

Figure  3  depicts  the  tantalum  tube,  heat  shield,  and  associated  plumbing  used  as  a 
variable  temperature  beam  source.  A  resistance-heated  tantalum  tube,  8  in.  long  and  3/16 
in.  in  diameter  with  a  wall  thickness  of  0.01  in.,  was  used  to  ensure  high  gas  stagnation 
temperatures.  The  center  of  the  tube  contains  a  0.01 4-in. -diam  hole  orifice.  Additional 
metal  was  removed  in  the  neighborhood  of  the  orifice  hole  to  ensure  a  thin-wall  orifice 
nozzle.  Large  copper  electrodes  attached  near  the  tube  ends  provided  electrical  power.  A 
Chromel®- Alum  el  ®  thermocouple  was  welded  directly  to  the  tantalum  in  the  vicinity  of 
the  orifice  and  was  used  to  monitor  source  temperatures.  A  water-cooled  copper  heat 
shield  surrounds  the  tantalum  tube  and  serves  to  protect  the  cryopanels  from  the  radiant 
heat  emitted  when  the  tantalum  is  used  as  a  high  temperature  source.  However,  liquid 
nitrogen  was  circulated  through  the  shield  when  a  low  temperature  source  was  desired. 
All  of  the  source  temperature  measurements  were  made  with  a  thermocouple  attached  to 
the  tube  behind  the  orifice. 

Z3  THE  SKIMMER  AND  COLLIMATOR 

A  conical  metal  skimmer  is  commonly  used  during  molecular  beam  experiments 
of  this  type.  However,  the  results  of  previous  investigations  (Ref.  10)  and  an  analysis  of 
the  flow  field  indicated  that  severe  skimmer  interactions  may  have  occurred  during  many 
of  the  previous  experiments  conducted  at  AEDC  and  other  laboratories.  Therefore,  the 
conventional  skimmer  was  removed,  and  a  copper  “doughnut”  cryocooled  by  gaseous 
helium  (22°K)  was  used  as  a  skimmer  during  the  experiments. 
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Figure  4  depicts  the  cooled  skimmer  used  during  these  tests.  The  skimmer  is 
constructed  of  1 /8-in. -thick  copper  sheet.  The  orifice  is  a  piece  of  l-in.-diam  copper 
tubing  1-3/4  in.  in  length.  The  outside  diameter  of  the  skimmer  was  tailored  to  fit  the 
hole  produced  by  removing  the  conventional  skimmer.  The  skimmer  outside  diameter  is 
approximately  8-1/2  in. 

At  low  flow  rates,  a  capture  coefficient  greater  than.  0.99  for  argon  on  a  22° K 
argon  cryofrost  is  reported  in  Ref.  11.  The  effects  of  high  flow  rates  and  changes. in 
cryofrost  crystallographic  structure  on  the  capture  coefficient  were  unknown.  An  RF 
antenna  was  installed  3  in.  in  front  of  the  skimmer  face  and  used  as  a  flow  field 
visualization  device  so  that  the  capture  ability  of  the  skimmer  and  degree  of  skimmer 
interaction  could  be  qualitatively  evaluated.  Figure  5  depicts  a  region  of  high  density  in 
the  flow  field  around  the  skimmer  at  room  temperature  (no  cryocooling).  Once  the  flow 
of  gaseous  helium  was  established  and  an  initial  argon  cryofrost  developed,  this  high 
density  region  completely  disappeared. 

2.4  THE  DETECTOR  SYSTEMS 

The  measurement  of  relative  abundances  of  molecular  clusters  and  the 
measurement  of  time-of-flight  (TOF)  data  require  different  types  of  detector  systems. 
Each  detector  configuration  is  discussed  below. 

2.4.1  The  Total  Flux  Detector  System 

Total  incident  beam  flux  was  measured  with  a  miniaturized  ionization  gage 
mounted  on  a  probe.  This  gage  could  be  accurately  positioned  in  the  flow  field.  The  gage 
was  enclosed  in  the  standard  glass  envelope,  but  the  normal  l-in.-diam  opening  was 
reduced  to  6  mm  as  shown  in  Fig.  6.  Reducing  the  opening  served  to  increase  the 
sensitivity  of  the  gage  to  directed  gas  flows  and  also  ensured  that  molecular  clusters  were 
broken  up  by  wall  collisions  so  that  a  gage  response  could  be  obtained  before  the 
molecules  escaped. 

The  total  beam  flux  was  obtained  from 


N  -  ai  (Pg  -  Pbg) 


where  Pb  is  the  measured  beam  pressure  and  Pbg  is  the  measured  background  pressure. 
The  calibration  constant,  «i  (2.5  x  1020  molecules/sec  torr  for  argon),  was  obtained  by 
using  an  oven-type  molecular  beam  to  provide  a  known  beam  flux. 

24.2  The  Modulated  Beam  Detector 

The  modulated  beam  detection  system  used  consisted  of  (1)  a  mechanical  beanr 
chopper,  (2)  a  quadrupole  mass  spectrometer  (2  to  600  amu),  and  (3)  a  lock-in  amplifier 
Figure  7  indicates  schematically  the  interconnection  of  these  components. 
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The  collimated  molecular  beam  is  chopped  into  pulses  by  two  slots  in  a  rotating 
wheel.  These  slots  are  90-deg  sectors  cut  out  of  a  2-in.-diam  disc.  The  beam  chopper  is 
driven  by  a  small  two-phase  synchronous  vacuum  rated  motor  generally  turning  at  40  Hz. 
These  beam  pulses  then  “fly-thru”  the  ionization  region  of  a  mass  spectrometer,  and  an 
ionized  sample  of  the  beam  is  electronically  deflected  through  the  quadrupole  section 
which  is  tuned  to  the  mass  number  of  the  beam  molecules.  The  resulting  ion  current  is 
amplified  by  an  electron  multiplier  which  sends  a  pulsed  electrical  signal  to  the  “lock-in” 
amplifier.  This'  amplifier  operates  as  a  bandpass  amplifier,  centered  on  the  chopper 
frequency.  A  light  and  photocell  mounted  at  the  chopper  wheel  provide  the  chopping 
frequency  to  the  amplifier  and  keep  it  synchronized  with  the  actual  molecular  beam 
modulation  frequency.  The'  lock-in  amplifier  increases  the  signal-to-noise  ratio  by 
amplifying  only  the  proper  frequency  and  phase  with  respect  to  the  reference  signal  from 
the  photocell.  The  amplifier  signal  is  displayed  on  an  oscilloscope,  on  a  chart  recorder, 
and  on  a  voltmeter  which  is  an  internal  part  of  the  amplifier.  Very  detailed  explanations 
of  this  particular  detection  system  are  contained  in  Refs.  9  and  12. 

Z4.3  The  Velocity  Distribution  Detection  System 

Reference  13  provides  a  detailed  description  of  the  detection  system  used  to 
measure  velocity  distributions.  A  summary  of  the  material  developed  in  the  cited 
reference  follows. 

The  detection  system  was  configured  as  illustrated  in  Fig.  7.  The  collimated 
molecular  beam  passes  through  the  narrow  slits  of  a  chopper  wheel  which  has  been 
designed  to  provide  an  approximate  impulse  of  molecules  (Ref.  13). 

As  time  elapses,  the  molecules  spread  according  to  their  velocity  distribution.  The 
beam  then  passes  through  a  focusing  aperture  plate  which  reduces  the  beam  diameter  to  a 
size  smaller  than  that  of  the  smallest  ionizer  aperture.  The  spectrometer  is  operated  as  a 
flow  density  detector  and  measures  the  local  particle  density  in  the  ion  source  as  a  time 
dependent  function.  The  local  density  in  the  ion  source  is  directly  proportional  to  the 
arrival  rate  of  the  molecules  and  inversely  proportional  to  their  velocity. 

The  quadrupole  section  is  tuned  to  the  particular  mass  number  of  interest,  and 
the  resulting  ion  current  is  amplified  by  the  electron  multiplier.  Essential  to  the  operation 
of  this  system  is  a  Field  Effect  Transistor  (FET)  emitter  follower.  The  FET  emitter 
follower  acts  as  a  buffer  or  impedance  matching  device  between  the  required  high 
amplifier  input  impedance  at  the  electron  multiplier  output  and  the  low  impedance 
required  for  the  coaxial  cables  which  lead  out  of  the  chamber.  Signal  losses  are 
minimized  by  mounting  the  FET  emitter  follower  circuit  directly  on  the  quadrupole  tube 
inside  the  vacuum  chamber. 

The  spectrometer  signal  is  given  a  time  reference  by  means  of  the  photoelectric 
cell  arrangement  shown  in  Fig.  8.  A  synchronous  pulse  referenced  to  chopper  open  time 
is  added  to  the  spectrometer  signal  in  the  differential  preamplifier.  This  composite  signal 
represents  the  molecular  arrival  rate  which  is  time  dependent  referenced  to  a  particular 
"start  time”  which  is  the  time  the  chopper  shutter  opens  plus  the  random  noise  of  the 
background  gas.  The  signal  is  then  processed  in  the  waveform  eductor,  a  schematic  of 
which  is  shown  in  Fig.  9. 
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Figure  9  traces  the  path  of  the  incoming  signal  through  the  waveform  eductor. 
The  incoming  composite  signal  is  fed  through  one  of  several  resistors  onto  a  signal  bus 
common  to  a  memory  bank  of  100  capacitors.  These  capacitors  are  connected  to  the 
signal  bus  by  means  of  an  insulated  gate-type  FET.  A  precision  oscillator  acts  as  a  clock 
and  advances  a  ring  counter  which  controls  the  memory  gates.  This  oscillator,  which  is  an 
internal  part  of  the  waveform  eductor,  is  externally  triggered  by  the  chopper  synchronous 
pulse.  The  effect  of  this  circuitry  is  to  allow  each  of  the  100  capacitors  to  be  connected 
for  consecutive  equal  time  periods  to  the  signal  bus.  By  synchronizing  this  switching 
action  with  the  chopper  and  the  flight  time  of  the  gated  molecular  beam,  the  same 
portion  of  the  velocity  distribution  is  applied  to  any  given  capacitor  on  successive  sweeps. 
Since  the  signal  is  being  integrated  as  a  result  of  the  selected  time  constant,  the  noise  is 
suppressed.  The  capacitors  are  also  connected  to  a  high  input  impedance  amplifier  which 
effectively  isolates  them  from  the  output  terminations.  The  memory  can  be  erased  or 
read  out  at  a  speed  which  is  compatible  with  the  data  recording  equipment  used.  The 
readout  speed  is  independent  of  the  speed  used  to  store  data. 

The  waveform  eductor  was  connected  to  three  output  devices: 

1.  An  oscilloscope  which  was  used  to  monitor  data  acquisition, 

2.  An  X-Y  plotter  which  served  to  maintain  a  permanent  record  of  the 
experiments,  and 

3.  A  data  logger  that  provided  a  punched  tape  which  was  used  as  a 
computer  input  for  data  reduction  purposes. 

SECTION  III 

EXPERIMENTAL  PROCEDURES  AND  FACTORS  AFFECTING  THE  DATA 
3.1  INTRODUCTION 

Two  separate  sets  of  experiments  were  conducted.  Initially,  the  concentrations1 
of  the  monomer,  dimer,  and  trimer  molecular  clusters  were  determined  as  plots  of  source 
pressure  versus  signal  strength  for  source  temperatures  of  170,  300,  and  600° K.  Velocity 
distribution  measurements  were  made  at  various  source  pressures  and  temperatures  and 
the  flow  field  properties  (most  probable  velocity,  static  temperature,  speed  ratio,  etc.) 
were  determined. 

&2  SOURCE  CONDITIONS 

A  source  to  skimmer  entrance  distance  of  22.2  cm  was  used  for  all  of  the 
experiments.  Calculations  indicate  that  this  separation  distance  will  ensure  free  molecular 
flow  through  the  detector  over  the  range  of  experimental  conditions. 


Concentrations  refer  here  to  the  actual  mass  spectrometer  signal  obtained  for  each  polymer  at  various  source 
conditions. 
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All  values  of  po  were  measured  with  calibrated  aneroid-type  dial  gages.  From  0  to 
800  torr,  an  accuracy  of  ±1  torr  was  obtained  with  the  gage  used.  In  the  range  from  800 
to  5000  torr,  an  accuracy  of  ±10  torr  was  obtained  with  a  Kollsman  gage.  Above  5000 
torr,  accuracy  was  considerably  reduced,  and  the  margin  of  error  is  estimated  as  ±500 
torr  for  this  gage. 

A  Chromel-Alumel  thermocouple  was  welded  to  the  outer  wall  of  the  tantalum 
source  opposite  the  source  orifice.  Previous  experimental  results  had  shown  that  the 
temperature  of  the  outer  wall  of  the  thin  wall  (0.010  in.)  tantalum  tube  would  be 
essentially  the  same  as  the  stagnation  temperature,  T0 ,  of  the  gas  in  the  tube. 

The  experimental  conditions  and  procedures  used  in  each  of  the  two  sets  of 
experiments  are  developed  separately  below. 

a3  MEASUREMENT  OF  THE  CONCENTRATIONS  OF  MOLECULAR  CLUSTERS 

These  measurements  were  conducted  using  experimental  procedures  developed 
during  previous  AEDC  tests  (Ref.  14).  A  summary  of  these  procedures  is  as  follows. 

The  total  beam  flux  versus  p0  was  measured  using  the  miniaturized  ionization 
gage  discussed  in  Section  2.4.1. 

The  modulated  beam  detection  system  was  used  to  measure  the  concentrations  of 
the  various  size  molecular  clusters  at  source  stagnation  temperatures  of  170,  300,  and 
600°  K.  AH  external  spectrometer  controls  were  optimized  and  held  constant  during  the 
measurements.  The  mass  spectrometer  was  initially  tuned  to  argon  mass  40.  The 
spectrometer  signal  was  recorded  as  the  source  pressure  was  incrementally  increased. 
Successively  higher  mass  number  signals  were  recorded  by  tuning  the  spectrometer  to  the 
new  mass  number  and  incrementally  traversing  the  source  pressure  in  alternate  directions 
(i.e.,  one  set  of  data  was  taken  with  the  source  pressure  varying  from  high  pressure  to 
low  pressure  and  the  next  set  was  taken  with  source  pressure  varying  from  low  to  high). 
After  taking  the  measurements,  mass  40  was  remeasured  to  determine  if  the  mass 
spectrometer  sensitivity  had  changed  during  the  test. 

The  higher  mass  number  molecular  clusters  required  fairly  high  beam  intensities  in 
order  to  properly  phase  the  lock-in  amplifiers.  Lock-in  was  accomplished  by: 

1 .  Raising  the  source  pressure  until  a  visible  signal  was  detectable  on  the 
oscilloscope, 

2.  Adjusting  the  detector  system  to  respond  to  the  species  of  interest, 
and 

3.  Incrementing  the  source  pressure  over  the  desired  range  and  recording 
the  mass  spectrometer  signal. 
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a4  FACTORS  AFFECTING  THE  DATA 

Interpretation  of  the  quadrupole  mass  spectrometer  signal  in  terms  other  than 
relative  signal  strength  involves  a  great  many  variables  which  are  not  all  well  understood 
at  this  time.  Since  the  purpose  of  the  molecular  cluster  concentration  measurements  was 
to  determine  the  source  pressures  and  temperatures  at  which  the  various  mass  number 
molecular  clusters  attained  a  maximum  value,  signal  strength  serves  as  an  adequate 
measurement  parameter. 

In  Refs.  3  and  14  are  discussed  inherent  errors  and  unknowns  in  the  use  of  mass 
spectrometers.  Those  areas  of  mass  spectrometry  which  influenced  the  experiment  are 
discussed  below. 

The  ionization  cross  section  of  argon  molecular  clusters  is  unknown.  Previous 
investigators  assumed  that  the  field  of  interaction  for  an  electron  is  small  with  respect  to 
the  diameter  of  an  argon  atom  and  that  an  electron  would  have  an  ionization  probability 
for  a  molecular  cluster  equal  to  the  sum  of  the  ionization  probabilities  for  the  atoms  in 
the  duster  (Ref.  14). 

The  possible  effect  of  cluster  fragmentation  due  to  the  ionization  process  in  the 
mass  spectrometer  is  also  difficult  to  assess.  During  these  experiments,  the  ionizing 
electron  energy  was  varied  from  20  to  90  ev,  and  the  relative  concentrations  of  the 
molecular  clusters  did  not  change  with  electron  energies  less  than  60  ev. 

The  electron  multiplier  sensitivity  is  a  function  of  the  age  and  environment  of  the 
multiplier  unit.  A  new  electron  multiplier  has  maximum  sensitivity  during  about  the  first 
eight  hours  of  use.  This  initial  eight  hour  period  is  also  the  time  when  its  sensitivity 
decreases  at  the  greatest  rate.  A  new  electron  multiplier  was  installed  in  the  mass 
spectrometer  at  the  start  of  the  molecular  cluster  concentration  measurements.  The  new 
multiplier  ensured  optimum  sensitivity  at  the  expense  of  repeatability.  A  recheck  of  data 
points  at  the  conclusion  of  each  experiment  produced  good  repeatability  but  overall 
spectrometer  sensitivity  decreased  by  approximately  1 0  percent  during  the  course  of  the 
experiment. 

In  summary,  the  goal  of  this  series  of  experiments  was  to  determine  the  source 
pressures  at  which  the  various  size  molecular  clusters  achieved  maximum  signal  values. 
With  this  objective  in  mind,  it  seems  reasonable  to  measure  concentrations  in  terms  of 
uncorrected  spectrometer  readings. 

3.5  MEASUREMENT  OF  THE  VELOCITY  DISTRIBUTION  OF  MOLECULAR 
CLUSTERS 

The  velocity  distribution  detection  system  was  installed  as  shown  in  Fig.  8,  and 
the  pressure  and  temperature  of  the  source  were  controlled  in  the  same  manner  as  before. 

Initial  phasing  of  the  detection  system  was  accomplished  by  operating  at  a  fairly 
high  source  pressure  and  tuning  the  quadrupole  section  of  the  spectrometer  to  the  desired 
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molecular  cluster  mass  number.  This  high  intensity  signal  could  be  monitored  on  the 
oscilloscope  and  readily  identified  by  its  chopped  appearance. 

After  phasing  was  completed,  the  source  pressure  was  reduced  and  TOF 
measurements  were  made  over  the  range  of  source  pressures  and  temperatures  by  the 
following  procedure: 

1 .  The  source  pressure  and  temperature  were  stabilized  at  desired  values 
and  the  miniaturized  ionization  gage  was  raised  out  of  the  path  of  the 
molecular  beam. 

2.  The  waveform  eductor  was  operated  in  analyze  mode  with  a  20-sec 
time  constant. 

3.  The  signal  storage  of  the  waveform  eductor  was  monitored  on  an 
oscilloscope  and  sufficient  time  (30  to  40  sec)  was  allowed  to  ensure 
complete  signal  storage  in  the  eductor.  This  complete  signal  storage 
was  attained  when  no  further  signal  increase  could  be  detected  with 
the  oscilloscope  operated  in  an  expanded  range. 

4.  The  eductor  was  then  operated  in  readout  mode,  and  an  X-Y  plot  and 
a  data  log  tape  of  the  stored  signal  were  obtained. 

5.  The  miniaturized  ionization  gage  (Fig.  6)  was  lowered  into  the  path 
of  the  molecular  beam  thus  turning  the  beam  “off.” 

6.  The  background -signal  was  stored  in  the  eductor  and  then  read  out 
on  the  X-Y  plotter  and  the  data  tape.  This  background  signal 
(beam-off  signal)  established  trie  base  line  reference  for  each  eductor 
module. 

Details  pertinent  to  TOF  analysis  are  included  in  Ref.  13.  The  computer  program 
used  to  process  the  data  tape  is  discussed  in  Ref.  1 5. 

The  functions  of  the  signal  distribution,  Sd(t),  computed  from  the  data  are 
affected  by  four  factors: 

1.  Accuracy  of  the  flight  length,  L,  measurement, 

2.  Accuracy  of  the  time  measurements, 

3.  Gas  purity,  and 

4.  The  agreement  of  the  fit  between  the  experimental  signal  distribution 
and  the  analytical  expression : 


Sd(t)  =  Ct~4  exp£-/S2  L2  - j 


(12) 
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which  is  derived  in  Ref.  15. 

Each  of  the  four  factors  is  discussed  in  turn  below. 

3L5.1  Accuracy  of  the  Flight  Length,  L,  Measurement 

To  minimize  the  effects  of  flight  distance  error,  it  is  desirable  that  the  flight 
distance  be  as  long  as  possible  within  the  physical  constraints  of  the  molecular  beam 
chamber.  Initially,  a  flight  distance,  L,  of  30.2  cm  was  used.  During  later  experiments, 
this  distance  was  increased  to  33.5  cm  by  an  improved  chopper  mounting  technique.  The 
accuracy  of  the  flight  distance  measurement  is  estimated  to  be  ±0.05  cm. 

3.1x2  Accuracy  of  the  Time  Measurements 

The  time  measurements  are  each  subject  to  different  degrees  of  error: 

1.  The  time  delay  between  the  photoelectric  trigger  and  the  pulse  of 
molecules  leaving  the  chopper,  tc,  was  obtained  by  comparing  the 
light  pulse  and  chopped  beam  signals  on  separate  traces  of  the 
oscilloscope.  This  measurement  was  made  to  an  accuracy  of  one 
microsecond,  and  its  error  is  considered  negligible. 

2.  The  time  delay  between  ionization  and  signal  detection,  tSl  was 
measured  by  using  a  precision  5-^ sec  pulse  generator  to  gate  the  ion 
beam  at  the  focus  aperture  of  the  mass  spectrometer.  The  5 -usee 
pulse  and  the  detected  signal  were  compared  on  separate  traces  of  the 
oscilloscope,  and  this  time  delay  error  is  considered  negligible. 

3.  The  time  delay  between  the  photoelectric  trigger  and  the  “turn-on” 
of  the  first  eductor  module,  to,  is  controlled  by  a  precision 
potentiometer  which  is  calibrated  in  1-usec  increments.  Calibration  of 
this  control  against  a  crystal  reference  standard  showed  a  maximum 
error  of  2  usee  when  time  delays  in  the  500-usec  range  were  used.  A 
500-usec  time  delay  is  a  typical  time  delay  used  in  the  experiments, 
and  this  error  is  considered  negligible. 

4.  The  sweep  duration  of  the  waveform  eductor,  100  At,  has  the  least 
precision  of  any  of  the  time  factors  discussed.  The  waveform  eductor 
sweep  duration  dial  setting  is  calibrated  in  50-usec  increments,  and  it 
is  possible  to  misalign  the  dial  by  ±5  usee  even  if  it  were  in  perfect  cali¬ 
bration,  When  calibration  error  is  included,  this  time  can  be  in  error  by 
as  much  as  ±  1 0  usee. 

3.5.3  Gas  Purity 

All  bottles  were  certified  as  at  least  99.995-percent  pure.  However,  one  set  of 
data  was  discarded  because  chemical  analysis  of  the  contents  of  one  gas  bottle  revealed 
gas  impurities.  After  this  experience,  the  contents  of  each  bottle  of  argon  were  analyzed 
by  traversing  the  signal  spectrum  of  the  mass  spectrometer. 
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3.5.4  Agreement  between  the  Experimental  and  Theoretical  TOF  Distributions 

The  single  factor  which  has  an  overwhelming  effect  on  the  data  is  the  quality  of 
fit  between  the  analytical  function  and  the  observed  data.  All  of  the  values  determined 
by  TOF  analysis  are  functions  of  the  analytical  expression,  Eq.  (12),  which  is  fitted  by 
the  least-squares  technique  to  the  experimental  data.  The  velocity  at  which  maximum 
signal  occurs,  vm ,  can  be  determined  with  a  high  degree  of  repeatability  since  by 
definition: 

1,  Flight  Distance 

m  tm  Time  of  Maximum  Signal 


The  parameter,  vm ,  is  least  influenced  by  the  quality  of  the  fit  since  it  is  dependent  only 
on  those  data  points  in  the  immediate  vicinity  of  the  point  of  maximum  signal  (Ref.  1 5). 
The  speed  ratio,  SR,  on  the  other  hand,  depends  on  ail  the  data  points  in  the  distribution 
curve  (Appendix  II)  and  is  thus  more  sensitive  to  variations  in  the  signal. 

SECTION  IV 

THE  EXPERIMENTAL  RESULTS 

The  results  presented  in  this  section  are  representative  of  molecular  cluster 
concentration  data  and  of  more  than  350  TOF  measurements.  The  polymer  concentration 
data  are  discussed  first  and  then  the  properties  (most  probable  velocity,  static 
temperature,  and  speed  ratio)  of  the  free-jet  flow  field. 

4.1  CLUSTER  CONCENTRATION  DATA 

Figures  10,  11,  and  12  denote  the  relative  signal  strength  produced  by  the 
molecular  clusters  as  a  function  of  source  pressure  (10  to  10,000  torr)  for  various  source 
temperatures  (170,  300,  and  600°K).  The  total  beam  flux  in  each  case  was  fitted  to  the 
mass  40  concentration  curve  in  the  region  where  no  dimers  were  detected. 

In  a  free-jet  expansion  of  argon  through  a  small  nozzle  (0-014-in.-diam),  three 
basic  processes  occur  as  the  source  pressure  is  varied  from  10  to  10,000  torr  for  a 
constant  source  temperature.  First  a  rapid  increase  in  the  dimer  concentration  is 
observed.  The  local  maximum  in  the  dimer  concentration  curve  occurs  at  successively 
higher  pressures  (90,  1500,  and  6000  torr)  as  the  source  temperature  is  increased  (170, 
300,  and  600°  K,  respectively)  as  shown  in  Figs.  10,  11,  and  12.  Second,  stable  low  mass 
unit  polymers  (e.g.  120  AMU)  appear  and  then  concentration  increases  very  rapidly  with 
increasing  source  pressure.  The  local  maxima  in  the  polymer  peaks  also  occur  at  higher 
pressures  as  the  source  temperature  is  increased.  The  low  mass  unit  polymers  decrease  as 
classical  condensation  begins.  However,  in  Figs.  10  and  11,  the  mass  80  and  120 
concentration  curves  begin  to  increase  for  increasing  source  pressures  after  the  onset  of 
classical  condensation. 
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42  FLOW  FIELD  PROPERTIES 

The  flow  field  properties,  i.e.  the  velocity  corresponding  to  the  time  of  maximum 
signal,  vm,  the  static  temperature,  Ts,  and  speed  ratio,  SR,  obtained  from  the  TOF 
measurements  are  presented  in  Figs.  13,  14,  and  15  as  a  function  of  source  pressure  for 
various  stagnation  temperatures. 

42.1  The  Velocity  Corresponding  to  the  Time  of  Maximum  Signal 

The  velocity  corresponding  to  the  time  of  maximum  signal,  vm ,  is  the  most 
accurately  determined  experimental  flow  field  property  and  is  defined  as 

L 

vm  =  7 
tm 

where  L  is  the  flight  distance  and  tm  is  the  time  of  the  TOF  maximum  signal.  For 
aerodynamic  molecular  beams,  i.e.,  high-speed  ratio  beams,  vm  is  approximately  the  mean 
velocity  of  the  flow  field  (Ref.  15).  In  Fig.  13a,  vm  for  both  mass  40  and  80  continually 
increases  as  P0  increases  for  a  stagnation  temperature  of  170°K.  However,  as  shown  in 
Fig.  14a,  vm  for  the  mass  40  signal  at  T0  =  300°K  increases  abruptly  at  a  source  pressure 
of  about  900  torr  and  then  shows  only  a  very  slight  increase  with  increasing  source 
pressure  beyond  1800  torr.  The  vm  for  the  mass  80  and  120  clusters  is  less  than  that  of 
mass  40,  but  a  slight  increase  with  increasing  p0  was  observed. 

The  energy  equation  for  a  perfect  gas  expanded  as  a  free  jet  to  very  small  static 
temperatures  is  closely  approximated  by 

cp  To  =  Y  m  V  2 

For  a  source  temperature  of  300°  K,  the  stream  velocity  for  a  free  jet  of  argon  is  5.5  x 
104  cm/sec.  As  shown  in  Fig.  14a,  the  measured  vm  is  approximately  5.45  x  104  cm/sec, 
and  when  condensation  occurs,  the  value  increases  to  approximately  6.2  x  104  cm/sec. 
On  a  macroscopic  scale  the  increase  in  the  mass  40  vm  at  approximately  900  tore  is 
possibly  the  result  of  the  latent  heat  of  condensation.  Based  on  this  hypothesis,  one 
theorizes  that  the  static  temperature  would  increase  (Fig.  14b)  when  the  latent  heat  of 
condensation  is  added  to  the  flow  field.  Thus,  more  molecular  collisions  would  occur  in 
the  expansion.  The  result  is  that  the  lower  molecular  weight  particles  (i.e.,  mass  40)  are 
accelerated.  The  exact  collision,  condensation,  and  evaporation  processes  that  occur  are 
not  completely  understood  but  the  preceding  provides  a  hypothetical  macroscopic  model 
of  the  general  processes. 

Figure  15a  presents  the  velocities  for  mass  40,  80,  and  120  as  a  function  of 
source  pressure  at  T0  =  600°K.  The  velocities  for  all  mass  numbers  were  observed  to 
decrease  as  p0  was  increased  from  400  to  6000  torr.  However,  beyond  6000  tore  (the  p0 
corresponding  to  the  onset  of  massive  condensation  in  Fig.  1 2),  vm  for  mass  40  began  to 
level  out  and  finally  slightly  increase  at  10,000  tore.  Unfortunately,  the  vacuum  system 
would  not  allow  increasing  p0  beyond  1 0,000  tore,  and  thus  the  existence  of  the  rise 
observed  in  vm  at  lower  temperatures  was  not  able  to  be  ascertained. 
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4.2.2  The  Static  Temperature 

The  static  temperature  of  the  flow  field  is  determined  by  using  a  least-square 
technique  (Ref.  15)  to  curve  fit  the  data  to  a  modified  Maxwellian  velocity  distribution 
function,  and  then  the  corresponding  temperature  is  calculated.  Representative  Ts  data 
for  source  temperatures  of  170,  300,  and  600° K  are  presented  in  Figs.  13b,  14b,  and 
15b.  For  all  cases,  the  static  temperature  at  first  decreases  with  increasing  source  pressure 
to  a  minimum  (below  1°K  for  T0  =  170  and  300°K).  Then  appreciable  condensation 
begins,  and  the  heat  of  condensation  is  released  into  the  flow  field.  The  result  is  a  rapid 
increase  in  the  static  temperature  with  increasing  source  pressures. 

The  static  temperature  data  are  the  least  accurate  data  presented  because  a  small 
error  in  curve  fitting  can  result  in  an  appreciable  error  in  the  value  of  Tg.  However,  the 
general  trends  are  repeatable  over  all  source  conditions. 

4.2.3  The  Speed  Ratio 

The  speed  ratio,  SR,  is  defined  as  U0A/2RTS  and  was  obtained  from  the  velocity 
distribution  data  as  discussed  in  Ref.  1 5.  The  data  for  various  values  of  T0  are  shown  in 
Figs.  13c,  14c,  and  15c.  All  of  the  data  exhibit  the  same  trends,  i.e.,  at  first  (as  pe 
increases)  there  is  an  increase  in  SR  to  a  maximum  value  and  then  there  is  an  abrupt 
decrease  with -increasing  source  pressure.  This  decrease  could  be  due  to  the  higher  static 
temperatures  that  seem  to  be  associated  with  the  onset  of  classical  condensation  (Figs. 
13b,  14b,  and  15b). 

The  maximum  speed  ratio  of  27.5  obtained  for  a  T0  =  300°K  and  p0  =  400  torr 
(Fig.  1 4c)  is  greater  than  the  maximum  predicted  by  Scott,  et  al.  (Ref.  1 6),  and  Hagena, 
et  al.  (Ref  17).  Condensation  does  decrease  the  speed  ratio  as  predicted  (Ref.  16),  but 
the  decrease  occurs  only  where  appreciable  condensation  is  present.  Skimmer  interaction 
has  been  virtually  eliminated  by  cryogenic  cooling  in  these  tests.  This  experimental 
condition  perhaps  accounts  for  the  fact  that  the  measured  speed  ratios  are  higher  than 
those  observed  during  previous  experiments  by  other  investigators  (Refs.  1 6  and  1 7). 

SECTION  V 
CONCLUSIONS 

This  series  of  tests  has  established  the  following: 

1.  Condensation  does  occur  in  free-jet  expansions,  and  it  does  affect  the 
properties  of  the  flow  field. 

2.  Condensation  tends  to  increase  the  mean  velocities  of  the  mass  40 
species  free-jet  expansions. 

3.  The  static  temperature  increases  when  appreciable  condensation 
occurs. 
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4.  Static  temperatures  lower  than  1°K  can  be  achieved  in  free-jet 
expansions  of  argon. 

5.  When  massive  condensation  occurs,  the  speed  ratio  is  greatly  reduced. 
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I.  ILLUSTRATIONS 
II.  THE  SPEED  RATIO 
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Fig.  3  Photograph  of  the  High  Temperature  Beam  Source 
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Fig.  4 


22  K  Skimmer 


of  the  Gaseous-Helium-Cooled  Skimmer 


p  =  725  torr  =  1270  torr 

°  n  ° 

T  -  300  K 
o 

Fig.  5  Flow  Field  around  the  Uncooled  Skimmer 
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Fig.  6  Total  Incident  Beam  Flux  Gage 
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Fig.  7  Schematic  Diagram  of  the  Modulated  Beam  Detection  System 
Used  in  Relative  Cluster  Abundances  Measurements 
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Fig.  8  Schematic  Diagram  of  the  Modulated  Beam  Detection  Systems 
in  the  Time*of- Flight  Measurements 
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Fig.  9  Schematic  Diagram  of  the  Waveform  Eductor 
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Fig.  10  Argon  Cluster  Concentrations,  T( 
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Fig.  11  Argon  Cluster  Concentrations.  T 
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Fig.  12  Argon  Cluster  Concentrations,  T0  =  600°  K 
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Speed  Ratio,  SR 


a.  vm  versus  p0 


c.  SR  versus  p0 

Fig.  13  Flow  Field  Properties,  T0  =  170°  K 


c.  SR  versus  Po 

Fig.  15  Flow  Field  Properties,  Tc  =  600°  I 
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APPENDIX  II 
THE  SPEED  RATIO 


The  speed  ratio,  SR,  is  a  parameter  used  as  a  descriptive  property  of  molecular 
beams  in  rarefied  gas  dynamics  literature.  It  is  defined  as 

SR  -  vSfr  tIM) 


In  Eq.  (II- 1),  U0  is  the  flow  velocity  of  the  beam,  T  is  the  static  temperature,  and  R  is 
the  gas  constant  per  unit  mass. 


Consider  the  denominator  of  the  speed  ratio 

v/2RTT  =  =  vmp  (II-2) 

where  vmp  is  the  most  probable  velocity  of  a  Maxwellian  gas  at  temperature,  Ts.  The 
speed  ratio,  SR,  can  be  thought  of  as  a  measure  of  how  well  the  expansion  process 
converts  the  random  thermal  velocity  of  the  gas  in  the  source  into  a  flow  velocity,  U0 . 


At  the  very  low  static  temperatures  associated  with  the  free-jet  expansion,  the 
energy  equation  for  a  perfect  gas 


cpTa  —  CpTs  +  —  mU0 

(H-3) 

closely  approximates  the  condition 

cpT0  »  -jpni  Uc2 

(II-4) 

and 

UQ  =  v'2cpT0/m  -  y}[2 y/(y  -  1)]  RT0 

(II-S) 

Normally,  an  increase  in  stagnation  temperature,  T0,  should  always  produce  an 
increase  in  the  speed  ratio,  SR.  This  relationship  is  drastically  altered  in  a  free-jet 
expansion.  A  few  nozzle  diameters  downstream  from  the  orifice  the  flow  becomes  free 
molecular,  resulting  in  a  “freezing”  of  both  the  Mach  number  and  the  temperature.  This 
highest  attainable  Mach  number  is  designated  the  terminal  Mach  number,  Mj,  and  it  is 
related  to  the  speed  ratio  by  Eq.  (11-6): 


SR 


u.  =  JT_  Uc 

v^nr;  vT~  vVhtb 


Vy72  (— )  MT 

V  a  /  V  2 


(II-6) 


Note  that 


SR  =  SR  (y,  MT) 

It  is  desired  to  measure  the  behavior  of  the  various  size  molecular  clusters 
throughout  the  flow  field  in  terms  of  a  single  parameter.  The  values  of  7  appropriate  for 
molecular  clusters  are  unknown.  If  values  for  the  clusters  could  be  determined,  7  could 
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still  not  be  treated  as  a  constant.  The  value  of  7  will  change  in  the  transition  region 
between  continuum  and  free  molecular  flow.  The  speed  ratio,  SR,  provides  a  means  of 
dealing  with  changes  in  M7  and  7  for  all  size  clusters  throughout  the  flow  field  in  terms 
of  a  single  concept. 


Additionally,  Anderson,  et  al.  (Ref  18),  reasoned  that  the  expansion  process, 
which  is  governed  by  the  stagnation  conditions  and  the  sonic  orifice  geometry,  should  be 
a  function  of  the  source  Knudsen  number,  Kns.  They  were  able  to  experimentally  verify 
that  argon  gas  follows  the  empirical  relationship 

i-y 

V1T  =  Kns  y  01-7) 


where 


(II-8) 


and  Dn  is  the  orifice  diameter  and  XD  is  the  mean-free  path  of  the  gas  in  the  source  at 
stagnation  conditions.  This  relationship  is  valid  over  the  entire  range  of  source  conditions 
reported  in  this  study. 


Empirically  then,  the  speed  ratio,  can  be  linked  to  source  conditions  by  the 
following  development: 


SR 


A 


SR 

SR 


-V? 

l-  Y 

Mt  y 

(II-9) 

1 

1  1  kT„ 

(II-10) 

=  vT 

4nr  n  4ur  Po 

-V? 

1  -  y 

(  l  1  kT 0  1  Y 

\  Dn  V 2  4rrr‘Po  J 

(11-11) 

1 

«  — - - 

1  T  V 

jj-j-  (constant  p0) 

(III  2) 

and  SR  will  decrease  as  T0  increases. 

To  summarize  the  development  thus  far,  it  can  be  said  that  the  speed  ratio  is  a 
descriptive  property  of  a  molecular  beam  which  can  be  related  by  valid  empirical  data  to 
the  source  conditions  if  argon  is  used  as  the  test  gas. 
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